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^ . Abstract. We study the dynamics of a nonlinear one-dimensional disordered system 

I from a spectral point of view. The spectral entropy and the Lyapunov exponent are 

extracted from the short time dynamics, and shown to give a pertinent characterization 
of the different dynamical regimes. The chaotic and self-trapped regimes are governed 
' by log-normal laws whose origin is traced to the exponential shape of the eigenstates 

of the linear problem. These quantities satisfy scaling laws depending on the initial 
state and explain the system behaviour at longer times. 
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1. Introduction 

(N 

^ I The motion of non-interacting particles in disordered lattices has been intensely studied 



over the last decades. In one and two dimensions, and for a sufficient amount of disorder 
in three dimensions, it has been shown that the spreading of quantum wavepackets 
is suppressed, a phenomenon known as Anderson localization [H [2]. However, the 
celebrated Anderson model that leads to this prediction is a highly simplified model 
which in particular neglects particle-particle interactions, so that a crucial question is: 



^ I Do interactions destroy (or, on the contrary, enhance) Anderson localization? 



A burst of interest on the subject was recently driven by the demonstration that 
these questions could be studied experimentally and theoretically with an unprecedented 
degree of cleanness and precision using ultracold atoms. This has produced an impressive 
number of results, concerning both Anderson localization (AL) [3], HI O |6l [TJ E] and the 
Anderson transition [9l HOl [HI [12], observed in 3D systems. Moreover, systems of 
ultracold bosons turned out to be very well modeled by mean-field approaches [T3| [Ti] . 
in contrast to fermionic systems where such a simplification of the corresponding many- 
body problem is not possible. 

From the theoretical and numerical point of view, ultracold bosons in a ID optical 
disordered lattice can be described rather realistically by a simple generalization of the 
original Anderson model including a nonlinear term taking into account interactions in 
a mean-field description. Numerical studies of the model suggest that the long term 
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motion of the particles is subdiffusive [TJl UHl HZl HSl CHI EHl EI]- These studies pointed 
out in particular the central role of chaotic dynamics in the destruction of AL. 

This approach has however two important drawbacks. The first is that the timescale 
of subdiffusion is larger by several orders of magnitude than the timescale of the single- 
particle dynamics, implying very long computer calculations which make difficult a 
full study of the interplay between disorder and interactions. The second is that the 
nonlinearity leads to a strong dependence on the initial conditions which also makes 
it difficult to give a "global" characterization of the different dynamical regimes. In 
previous works, we have shown the existence of scaling laws with respect to the width 
of the initial state, allowing, to some extent, such a global characterization j2l], and 
demonstrated that these scaling laws are robust with respect to decoherence effects |22] 
(which can also destroy AL). 

In the present work we tackle the problem of interacting ultracold bosons in a ID 
disordered lattice using a spectral analysis that does not require such long computational 
times, as the information on the chaotic behaviour is "inscribed" even in early times in 
the spectrum of the dynamics. This allows us to perform a more complete and precise 
study of the problem over a large range of parameters. A central quantity in our 
study, the spectral entropy, proves very useful to characterize the dynamic behaviours, 
confirmed by comparing the information extracted from the spectral entropy to the 
Lyapunov exponent, another well-known measure of chaotic behaviour. Moreover, we 
show that such behaviours are described by "log-normal" laws which can be scaled with 
respect to the initial conditions, and we propose a simple physical interpretation of our 
findings. 

2. The model 

We use here the discrete nonlinear Schrodinger equation with diagonal disorder. It is 
essentially the same model used in previous works [2T1 [22], so we only give an outline 
of it here. 

The mean-field theory applied to ultracold bosons in an optical (ordered) lattice 
leads to the so-called Gross-Pitaevskii equation: 

(j){x) 

where is the macroscopic wavefunction (or order parameter) describing a Bose-Einstein 
condensate well below the critical temperature, V{x) = — Vq cos{2kLx) is the optical 
potential and qid the ID coupling constant [23j. Tight-binding equations are obtained by 
decomposing the wavefunction 0(x) = Yln'^n{t)wn{x) onto the set of localized Wannier 
functions of the first band w„(x) associated to the n^^ lattice site: 

VnCfi Cn—1 ~\~ Q \Cn\ Cj^, (l) 

where we kept only (symmetric) nearest-neighbours couplings, which is justified if the 
Wannier functions are strongly localized, that is, for large enough Vq. According to usual 



ih(f){x) 
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j- + V{x)+g^DN\<jy{x)\' 
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conventions, we measure distances in steps of the lattice and write energies in units of the 
couphng constant of neighbour sites T = — J dxWn{x)[p'^ /2m + V{x)]wn+i{x)^. Finally, 
times are written in units of h/T. The coefficient Vn = J dxWn{x)\p'^ /2m+V {x)]wn{x) /T 
is the diagonal on-site energy and the effect of interactions is taken into account, in 
a mean- field approach, by adding the nonlinear term g\cn\ where the dimensionless 
interacting strength is 



In the absence of disorder, the on-site energy f„ does not depend on the site index n and 
can be thus set to 0. According to the Anderson's postulate [IJ, we introduce diagonal 
disorder by picking random on-site energies f„ uniformly in an interval [—W/2, W/2]. 

For g = 0, describes the standard Anderson model, and we shall call the 
corresponding eigenstates (eigenvalues) "Anderson" eigenstates (eigenvalues). A few 
facts about it will be useful in what follows. The eigenstates are exponentially localized 
in average c'^ ~ exp {—\n\/l,y) (the overbar indicates the averaging over realizations of 
the disorder) with a locahzation length 1^{W) ~ 96(1 - e^/4)/VF^ [21]. For 51 7^ the 
equation becomes nonlinear, and it is useful, if not strictly rigorous, to interpret the 
nonlinear term as a "dynamical correction" to the on-site energy v^^ = glcnl"^- Previous 
works [T5| [T6| [TTf [20| [T8| 121] put into evidence the existence of three main dynamical 
regimes: For g <C W, Anderson localization is expected to survive for very long times, 
a regime that we shall call "quasi-localized". For g ~ W, the nonlinear correction 
v^^ induces chaotic dynamics leading to subdiffusion and the destruction of AL. For 
g ^ W, the very large nonlinear term v^^ decouples all sites whose populations are not 
nearly equal (even in the absence of the disorder), suppressing the diffusive behaviour 
and leading to another type of localization, called self-trapping. This dynamics does not 
rely on quantum interference and is therefore very robust against external perturbations, 
including decoherence |22) . 

Our aim is to characterize the global dynamics on a relatively short timescale. We 
thus study the evolution according to ([1]) of bosons in a ID box containing L sites 
(typically, L = 101) and put an exponential absorber at each end of the box in order 
to prevent wavepacket refiection. The norm of the wavepacket is thus not anymore 
conserved as soon as it "touches" the borders, and we characterize the diffusive behaviour 
by calculating the survival probability p{t) = |c„p; a value p(t) < 1 indicates that 
the packet has diffused outside the box. As in ref. [21], we restrict the analysis to initial 
wavepackets of the form 



§ As discussed in more detail in ref. [21 , this form of wavepacket has the advantage of, on the average, 
projecting onto all Anderson eigenstates, thus rendering the dynamics roughly independent of the 
wavepacket energy. 



9 = 



T 



with Lq <^ L 




I Wannier functions have the translation property Wn{x) = wq{x — n). 
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Figure 1. Survival probability p at time t = 10^ vs interaction strength g for W = 4 
and different widths of the initial state : Lq =3 (blue squares), 7 (green triangles), 13 
(red diamonds), 21 (cyan stars), 31 (magenta circles), 41 (yellow inverted triangles). 



The typical behaviour of the survival probability is illustrated in figure [T], which 
represents p{g,t = 10^) as a function of the interaction strength g for W = 4 and 
for various values of Lq. To obtain such smooth curves, the survival probability was 
averaged over typically 500 realizations of the disorder and of the initial phases 6'„. Such 
long times are necessary in order to clearly put into evidence the three different regimes 
mentioned above, which are then easily identifiable: At low g, the survival probability 
is close to one, corresponding to the quasi-localized regime in which AL survives. For 
intermediate values of g, the decrease of p indicates that the wave packet has spread 
along the box and part of it has been absorbed at the borders, corresponding to a 
chaotic regime induced by the nonlinearity that destroys AL. For large values of g, the 
survival probability increases again and gets close to 1, indicating that the packet is 
self-trapped [T5| [T6| 121] and has never touched the borders if initially it was thin enough 
(e.g. in the Lq = 3 case). In the next section, we show that a complementary way 
to describe the dynamics allows to a very precise characterization within much shorter 
computation times. 



3. Characterizing chaotic dynamics with the spectral entropy 



Spectral analysis is a very useful way of analyzing a chaotic dynamics, be it in classical, 
"quantum" [Jj] or quantum nonlinear chaotic systems [25]. The spectral entropy is a 
measure of the "richness" of a spectrum. Chaotic behaviours are associated to continuous 
spectra, and thus to a high spectral entropy. Given a quantity M{t), its power spectrum 
is defined as 
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II Traditionally the term "quantum chaos" designates the behaviour of a quantum (linear) system whose 
classical (nonlinear) counterpart is chaotic. 
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where M{f) is the Fourier transform of M{t) for t G [0,tjnax]- The choice of the value 
of ^max determines the resolution of the spectrum. In our case, we chose t^nax = 200, 
which is a very small value compared to the typical time of emergence of the interacting 
regime (t ~ 10^) but which will be shown to be sufficiently high to characterize the 
dynamics from a spectral point a view. As we do not expect excitations whose timescale 
is inferior to the tunneling time (1 in our rescaled units), we set /max = 10- From the 
power spectrum, one defines the spectral entropy as 

^ Jdf Sif)\ogS{f) 

log(/max) 

For a perfectly monochromatic signal Sf = 6{f — /o), the spectral entropy is zero, 
whereas for a white noise (5*/ = l//max) H = 1. The spectral entropy "counts" the 
number of frequencies which are present in the signal, and is a good indicator of the 
chaoticity of the system [26j. The spectral entropy obviously strongly depends on the 
choice of the observable, and its usefulness as a dynamics indicator is reliant on this 
choice. 

A good observable in the present problem is the so-called "participation number" 
P with respect to the Anderson eigenstates, which is defined, in the present case, as 
follows. For a given realization of disorder {f„}, we calculate the Anderson eigenstates 
in the Wannier basis, (f),^{x) = ^^^^dn^Wnix) corresponding to an energy e,^, which are 
solutions of ([T]) with g = (the d^^ replacing the c„). Back to the g case. Equation 
([1]) allows us to calculate the evolution of the wavepacket "ipit) = Cn{t)wn{x) under the 
action of both disorder and nonlinearitjQ. At any time, we can express the wavepacket 
in the Anderson eigenstates basis 

= $^?.W0.(a;), (3) 



from which, trivially, qu{t) = Yin d". c„(t). The participation number is then defined as: 

p = (4) 

U g = 0, the (f),y are the exact eigenstates of the problem, so that the populations \q^\'^ 
are constant. In the non-interacting case g ^ 0, the Iq^l"^ evolve under the action of 
the nonlinearity. The participation number roughly "counts" the number of Anderson 
eigenstates participating significantly in the dynamic^ and its time evolution thus 
refiects the apparition of Anderson eigenstates that were not initially populated. 

Figure [2] shows the spectral power S{f) of the participation number P in the three 
interacting regimes: (a) g = 1 in the quasi-localized regime, (b) g = 100 in the chaotic 
regime, and (c) g = 1000 in the self-trapping regime. For g = I, the dynamics is very 
similar to the linear case: The populations of Anderson eigenstates practically do not 

% None of our results is modified if we set /max = 2. 

+ We use a standard Crank-Nicholson scheme with time-step 0.01 < dt < 0.1. 

* To see this intuitively, consider two limit cases: If only one i/ = i/q Anderson eigenstate is populated, 
= 5u,uo ^"^"^ thus P = 1; if Lo eigenstates are equally populated, = ^ ^'^"^ ^ ~ ^o- 
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/ / / 

Figure 2. Example of spectral power S{f) for the three different regimes : a) localized 
g = 1, b) chaotic g — 100 and c) self-trapped g — 1000 . Other parameters are W — 3, 
Lq = 3 and f^ax = 200. 



evolve in time, as well as the participation number P and the spectrum is dominated 
by low frequencies. As a consequence, the spectral entropy, calculated according to (|2]), 
is relatively small: H = 10"^. For g = 100, Anderson eigenstates are strongly coupled 
by the nonlinear term, each pair of coupled states generating a Bohr frequency (shifted 
by the nonlinear correction), that is ~ Cj, + (7 ~ tyi — g \qui\^ ■ In this regime, most 
Anderson eigenstates are coupled to each other, so that the power spectrum is almost 
fiat (with important local fluctuations) at high frequencies, and the spectral entropy 
increases by almost an order of magnitude H = 10~^ with respect to the preceding 
case. For g = 1000 the wavepacket is self-trapped and only a relatively small number 
of Anderson eigenstates whose populations happen to be close enough can interact. 
The spectrum is therefore dominated by a finite number of frequencies and the spectral 
entropy is reduced, H = 2 x 10~^. However, although populations are stable in this 
case, quantum phases may evolve chaotically under the action of the nonlinearity |27| . 
Our definition of the spectral entropy from the participation number - which does not 
directly depends on the phases - excludes this "phase dynamics" from the corresponding 
spectrum. 

We display in figure [3^ the averaged spectral entropy as a function of the 
nonlinearity parameter g for different widths Lq of the initial state. One clearly sees the 
crossover from the quasi-localized to the chaotic regime, signalled by a marked increase 
of H. The smaller the value of Lq, the smaller the value of the crossover. This is easily 
understandable, as a more concentrated wavepacket leads to a stronger nonlinear term 
f^^. On the right side of the plot, one also sees, especially for low values of Lq the 
beginning of a decrease of H due to self-trapping. 

It is interesting to compare the information obtained from the spectral entropy to 
another relevant quantity characterizing chaos, the Lyapunov exponent, which indicates 
how exponentially fast neighbour trajectories diverge. This quantity is usually defined 
for classical systems, but can be extended, with a little care, for quantum nonlinear 
systems |[25j. We present in Appendix A a method for calculating the Lyapunov 
exponent of a quantum trajectory defined by amplitudes c„ [Equation ([1])]. Figure [3)d 
displays the Lyapunov exponent A obtained with the same parameters as in figure [3^.. It 
displays a monotonous increase with the nonlinear parameter, even in the region where 
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Figure 3. (a) Spectral entropy H and (b) Lyapunov exponent A vs interaction strength 
g ioi W — 4: and different widths of the initial state : Lq —3 (blue squares), 7 (green 
triangles), 13 (red diamonds), 21 (cyan stars), 31 (magenta circles), 41 (yellow inverted 
triangles) . 



the spectral entropy decreases due to self-trapping, evidencing the presence of a regime 
of "phase chaos" mentioned above. This clearly shows that H and A provide different 
information on the dynamics of the system, and that H is clearly more adapted to 
distinguish the three dynamical regimes. 

In section H] we will discuss the shape of these curves, their scaling properties, and 
suggest a physical mechanism explaining such properties. 



4. Log-normal law and scaling 



As we shall see below, log-normal functions are ubiquitous in the dynamics described 
in the present work, and so in various contexts which are not obviously related to 
each other. A first example can be found in figure [H where each curve representing 
p{logg) can be fitted with a rather good accuracy by a inverted Gaussian function. 
More generally, a log-normal function is defined as 

\2" 



fix) 
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exp 



(logo; - 
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exp(-/i + 
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exp 



(logX — /U + (T^)^"' 



In physics, such a function appears more often as a "log-normal distribution", related 
to the statistics of quantities which are a product of randomly distributed terms |28| . 
Log-normal statistics is very different from normal (Gaussian) statistics; for example, 
the most probable value of a log-normally distributed quantity is different from its 
average value. In the following, we shall not only consider statistical distributions over 
the realizations of the disorder: figure [1] does not display the distribution of p over the 
realizations of disorder but its average as a function of the interacting strength, so does 
figure [3^. for the spectral entropy. However and more interestingly, we now show that 
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Figure 4. Histograms of the spectral entropy H (a) and the Lyapunov exponent A 
(b) for = 3, Lq = 41 and two interacting strengths : 5 = 10 (left hand-side blue 
histogram), g = 100 (right hand-side red histogram). Black solid lines correspond to 
log- normal fits. 



the statistical distribution of the spectral entropy and of the Lyapunov exponent are 
indeed log- normal. 

Let us thus study the distribution of these two quantities over the realizations of 
the disorder u„ and of the initial phases 6n- Figure displays the distributions of 
values of H over 10^ realizations of the disorder for two values of g, and figure Hb the 
corresponding distribution for the Lyapunov exponent A. As shown by the black fitting 
lines, both curves are perfectly fitted by a log-normal function. 

In order to have an idea of the origin of these shapes, one can consider a simple, 
heuristic model. The destruction of AL is due to the nonlinear coupling between 
Anderson eigenstates. Initially unpopulated eigenstates, which, in the absence of 
nonlinearity, would never be populated, can be thus excited thanks to a nonlinear 
transfer of population. The goal of the simple model developed below is to characterize 
the statistical distribution of such excitations. Projecting the wavepacket in the 
Anderson eigenbasis [cf. Equation ([3])], ([T]) then reads: 

iqu = e^Qu + 9 ^li^^^ffyJi^^ ^2, z^a) (5) 

with /(z/, z/i, z/2, z/3) = dn^ dn^"* dn'^^ dn'^^ Populatious exchanges are controlled by the 
overlap / of four coefficients. The population transfer from eigenstate /i to eigenstate u 
depends on Ji = /(/i, /i, i/, z/), J2 = /(z^, t-') and J3 = /(/x, z/, z/, v) , the term due to 
Ji being for instance : 

^ = 2MIm(g:^gJ) 

In order to evaluate the probability distribution of Ji, we make the assumption that 
the coupled Anderson eigenstates are exponentially localized with the same localization 
length ^. The exponential localization is valid on the average and the localization length 
is the same if the eigenstates have close enough eigenenergies. We thus write 

4^)^ = tanh exp 
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= tanh (I) exp 




2\n — l{fi, z/)| 



) 




1 - e~4/« 



2 



+ /(/i,z/)-l 



(6) 



where, for simplicity, we have supposed that the spatial distance between the eigenstates, 
/(/X, u), is an integer. The most important term in ([6]) is e"^''^^''^-'/'^. In the limit ^ — )■ 0: 



where G = log [/(/i, u) + 1] — 2/(yU, z/)Ao and = 4:l{fi, vf'o'^ ■ The overlap sum J\ thus 
obeys a log-normal distribution. 

We conjecture that overlap sums like Ji, controlling the coupling between Anderson 
eigenstates, in fact control the destruction of the Anderson localization, and thus 
explain the log-normal shapes we observed for H and A. The above heuristic argument 
undoubtedly presents various assumptions that are not rigorously justified, but it has the 
merit of putting into evidence the intimate relation between the log-normal distribution 
and the exponential localization of the Anderson eigenstates. The link between the 
exponential shape and the emergence of log-normal statistics has also been studied in 
the case of the conductance of disordered systems |30[ [3T| |21] . 

Let us now consider the averaged spectral entropy. In previous works [2T1 |22] . 
we showed that suitable scaling with respect to the initial state width Lq allowed a 
classification of dynamic regimes independently of the shape of the initial state. In 
particular the interacting strength g was scaled as ^ = gL'^^ with s ~ 3/4. This scaling 
is meaningless for low values of Lq because in this case, the initial participation number 
is not of the order of Lq but is of the order of the maximum Anderson localization 
length l^iyV^ ~ 96/W'^. Figure [5^ shows that scaling with g and H = HL^ make the 
curves H{g) corresponding to Lq > 20 collapse to a single curve, except in the strong 
self-trapped region. Figure [Sjo shows the Lyapunov exponent as a function of g; the 
curves also collapse for Lq > 12. The Lyapunov exponent itself is independent of Lq 
(with no scaling of the variable A), which is not surprising as it does not measure 
the absolute distance between two quantum trajectories but the timescale of their 
exponential divergence. 



J^ = ^^2lM/^ [/(/i, z/) + l] 

The inverse localization length A = 1/^ follows a normal distribution [291 [21] 

P(A) oc exp[-(A - Ao)V2(t2] 
we obtain for the distribution of values of the overlap Ji 
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Figure 5. (a) rescaled spectral entropy H = HLq with s w 3/4 and (b) Lyapunov 
exponent A vs the rescaled interaction strength g ioi W = A and for Lq =3 (blue 
squares), 7 (green triangles), 13 (red diamonds), 21 (cyan stars), 31 (magenta circles), 
41 (yellow inverted triangles). 



As shown by the black sohd hne in figure [5^1, the scaled spectral entropy H is very 
well fitted (outside the strong self-trapped region g > 100) by the log-normal law: 



H 



hr 



gV2 



: exp 



2a2 



(7) 



with three free parameters, the amplitude Hq, the center G of the distribution and 
the "standard deviation" a. The study of these fitting parameters as a function of 
the disorder W provides a full characterization of the dynamic regime. Instead of 
representing the fit parameters and G, we prefer to use more physical quantities, 
namely the maximum value of H, -ffmax = ho (27ia'^) exp (o"^/2 — G) (figure [6^) 
and the rescaled interaction strength cjc = exp (G — cr^) (figure [Ub) corresponding to 
this maximum. The dependence of the standard deviation a on W is displayed in 
figure int. The quantity i^max is a decreasing function of W: as the localization length 
decreases with W, so does the overlap between two neighbour Anderson eigenstates. On 
the contrary, c/c is an increasing function of W: The number of resonances is maximum 
when v^^ is comparable to the typical energy between two neighbour states, itself of 
the order of the bandwidth ~ 4 + W; cjc is thus independent of W at low disorders and 
increases with W for larger disorders. Finally, in figure [Gt, one can notice that the log- 
normal curve becomes sharper when the disorder increases, which can be attributed to 
the fact that Anderson Localization is more robust against interactions at high disorders. 

The fact that the spectral entropy can be determined from a relatively short time 
interval also allows one to follow the evolution of the dynamics. Figure [7] shows the 
evolution of the spectral entropy H{t), defined as the spectral entropy calculated in an 
interval + 200] , that we shall call, for short, the dynamic spectral entropy. For 
the low value of (7 = 10 (figure [7^), after the destruction of AL, the diffusion of the 
wavepacket produces a dilution and a consequent diminution of the nonlinearity that 
reduces the chaoticity, and thus the spectral entropy of the system. For the high value 
oi g = 1000 (figure Wp) one sees a more complex interplay of different regimes: The 
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Figure 7. Evolution of the spectral entropy H{t) for W — 4, Lo = 41. For (a) g = 10 
one essentially sees the effect of dilution, which progressively weakens the chaoticity 
of the system. For (b) g = 10^, there is a first phase, in which the self-trapping 
is progressively destroyed which is followed by a slow transition towards the chaotic 
regime. 



initial state is initially frozen in a self-trapping regime, but this regime is unstable (for 
this particular set of parameters): If a fraction of the packet escapes the self-trapped 
region, the nonlinear contribution decreases which leads to a weakening of the 
trapping. The destruction of the self-trapping regime takes a much longer time that 
the destruction of the AL. One first observes a small decrease of the spectral entropy, 
as some eigenstates leave the box and do not interact anymore. Then, when v^^ has 
decreased sufficiently in the center of the packet, the system enters the chaotic regime. 

In plots a and b of figure [8] we represented the dynamic spectral entropy H{t) for 
three values of Lq, corresponding to a same value of g (0.62 in a, 6.17 in b and c). 
The curves converge for long times, putting into evidence the existence of a universal 
asymptotic regime, independent of the initial state, once the proper scaling on g is 
applied. Additional scaling can even be used to describe the transition from the self- 
trapping regime to the chaotic regime, as shown in figure IHfc. Although we have presently 
no physical explanation for the exponents appearing in these scaling laws, these results 
strongly support the idea that the features of the nonlinear dynamics should be scaled 
again with respect to the initial state [2T] . 
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Figure 8. Scaling of tlie spectral entropy H{t) ioT W — A and Lq = 13 (blue full 
line), Lq — 21 (green dotted) and Lq = 41 (red dashed). The use of the same values 
of the scaled nonlinearity g = gL^ '^ [g = 0.62 (a) and g = 6.17 (b)] shows that there is 
a universal asymptotic regime. Comparison of plots (b) and (c) shows that additional 
adhoc scaling gives an even more universal character to this behaviour. 



5. Conclusion 

We have shown that the spectral entropy is a very good indicator of the chaoticity 
of the dynamics, and allows a very good characterization of the dynamics regimes. 
Moreover, it can be calculated dynamically and also gives information on the evolution 
of the dynamics. The Lyapunov exponent gives, in the present context, a less complete 
characterization of the dynamics. It describes very well the progressive destruction of 
the Anderson Localization by the onset of chaotic behaviour, exactly as the spectral 
entropy, but it does not give information in the self trapping regime, where "phase 
chaos" is still present. 

The argument presented in section H] suggests that the log-normal shape is linked 
to the exponential localization of the Anderson eigenstates. It is a bit surprising to find 
this link between the behaviour of a strongly nonlinear system and the eigenstates of the 
corresponding linear system. From a mathematical viewpoint, spectral analysis (in the 
sense of the determination of eigenvalues and eigenvectors) cannot be directly applied to 
nonlinear systems, and, presently, there is no alternative analytical method for analyzing 
nonlinear systems. The results discussed in the present work and in refs. [2T| [22] indicate 
that scaling over the initial state is a promising tool for the study of the emerging field 
of nonlinear quantum mechanics. In this context, it is worth noting that we have 
presently no convincing explanation for the precise values of the exponents appearing 
in these scaling laws, and that we have considered only a particularly simple family of 
initial states. Moreover, there is no experimental verification of these scaling laws. It 
thus appears that a large field of investigations is opened for both theoreticians and 
experimentalists in the near future. 
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Appendix A. Lyapunov exponent of a quantum trajectory 

We are interested in calculating the Lyapunov exponent of the quantum trajectory 
c = (c„). Considering two initial conditions c", c*, which are separated by an 
infinitesimal distance do, the Lyapunov exponent measures the rate of their exponential 
divergence. 

A = lim- log ^ (A.l) 

t-i-co t do 



where d{t) = \c^{t) — c"(t)|. In principle, one could use directly (lA.ip but it turns out 



that the trajectories do not evolve in an ensemble of infinite volume. As a consequence 
the distance d{t) rapidly saturates and A tends to 0. The common method to get rid 
of this drawback is to let trajectories evolve during a short time period dt and then 
evaluate the corresponding Lyapunov exponent 

Ao = y log ^ (A.2) 

"0 

where di = \c^{dt) — c"'{dt)\. Before letting the system evolve for another time 
interval dt, we rescale trajectories so that the distance between a and b is set to do keeping 
their relative orientation unchanged. Usually, people rescale the second trajectory 
= + ^(c* — c^^jll satisfying immediately both requirements. We then obtain 
Ai, A2, .. by iterating the three-steps operation : (i) evolution during a short time dt (ii) 
calculation of the corresponding Lyapunov exponent from ( 1A.2I) (iii) renormalization of 
the trajectories. We finally deduce the final Lyapunov exponent : 

1 ^ 

A = lim — > Aj 

where A represents the number of iterations. Unfortunately, the crucial rescaling 
operation changes the norm of the second trajectory c;, : in the case of quantum systems, 
the procedure is therefore based on non-physical states. Here we propose to modify the 
step (iii) rescaling both trajectories using the following scheme : 

~e = ac" + /3c^ 
= + 6c^ 

where {a, (3, 7, 6) satisfy the following condition Jffl : 



d, ^ 
1 

1 



see for example 



http://sprott.physics.wisc.edu/chaos/lyapexp.htm 



(A.3) 

(A.4) 
(A.5) 



ft We calculate the Lyapunov exponent without using an absorber potential so that the norm |c| is a 
conserved quantity. 



Spectral description of the dynamics of interacting bosons in disordered lattices 14 



Given that di 



the first equation is nothing but the usual scahng condition 
which is used for classical systems. Projecting (1A.3P on and on c^, one immediately 
obtains 

do 
di 
do 



7 



a — 



5 = /3 + 



di 



Subtracting ( 1A.4I) from (1A.5I) . and using 23f?(c'^|c'') 
do 



2 — df, we obtain 



(3 



a — 



di 



and finally from (1A.4I) . we can choose 



a 




to finally deduce /3,7, 5. This method allows an accurate calculation of the Lyapunov 
exponent in the discrete system considered in the present work, but it can in principle 
be also applied to the continuous Gross-Pitaevskii equation. 
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